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FOREWORD 


The  investigation  described  in  this  report  was  conducted  by 
the  Aerodynamics  Branch  of  the  Aircraft  Laboratory,  Aeronautics 
Division,  Wright  Air  Development  Center.  The  work  was  initiated 
by  a Suborder  dated  12  December  1951  initiated  by  Mr.  J.  H.  Allen 
of  the  Equipment  Laboratory,  WADC.  The  CEO  No.  is  C-03. 

The  model  tests  were  performed  by  the  Dynamic  Model  Unit  under 
the  direction  of  Major  A.  J.  Stolzenberger.  The  performance  calcula- 
tions, evaluation  of  the  photographic  records  of  the  model  tests,  all 
the  drawings,  and  an  independent  check  of  all  numerical  calculations 
included  in  this  report  were  made  by  1st  Lt  E.  D.  Wong*  This  help 
is  gratefully  acknowledged  herewith  by  the  author. who  acted  as  project 
engineer. 
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ABSTRACT 


This  report  investigates  the  stability  of  a gondola  suspended 
from  the  "whirl  arm  of  the  Parachute  Test  Tower  at  El  Centro, 

California  (see  Fig,  l)  on  which  instability  had  been  observed* 

The  work  was  started  with  an  investigation  of  the  pendulum 
oscillations  of  the  gondola  neglecting  the  coupling  between  the  two 
modes*  .The  result,  that  both  frequencies  become  equal  at  approximately 
the  speed  at  which  the  instability  had  been  observed,  was  interpreted 
as  an  inclination  towards  instability*  Therefore  a more  detailed 
investigation  including  the  coupling  between  the  two  pendulum  modes 
was  made* 

The  detailed  investigations  showed  that  the  damping  of  the 
pendulum  modes  is  extremely  low.  Stability,  however,  is  secured  if 
the  gondola  has  longitudinal  static  and  weathercock  stability  and  if 
the  slope  of  its  cross  force  versus  angle  of  yaw,  , is  smaller 

than  the  drag  coefficient  of  the  gondola  plus  cable* 

Tests  made  on  a 1:20  scale  model  of  the  parachute  test  tower 
verified  the  essential  parts  of  the  theoretical  investigations,  namely 
the  predicted  dynamic-stability,  the  low  damping  and  the  order  of 
magnitude  of  the  frequencies* 
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SYMBOLS 


Aq,  A^,  A^,  j 

\ 0 

Coefficients  of  the  frequency  polynomial 

&L1> 

0 

Coefficients  of  the  constant  terms  in  the 
elements  of  the  frequency  polynomial 

bll»  b12> • • *b44 

0 

Coefficients  of  the  linear  terms  in  the 
frequency  determinant 

* 

(lb) 

Centrifugal  force 

C 

(lb) 

Aerodynamic  cross  force 

cll»  c12» • • *CA4 

0 

Coefficients  of  the  quadratic  term  in  the 
frequency  determinant 

c « c 

°c  /f  va® 

0 

Cross  force  coefficient  of  the  gondola 

Cd 

0 

Drag  coefficient  of  the  gondola  including 
part  of  the  cable  drag,  based  on  gondola 
cross-section. 

Cm 

0 

Pitching  moment  coefficient  of  gondola 

Cn 

0 

Yawing  moment  coefficient  of  gondola 

D 

(lb) 

Drag 

3 

(ft/sec2) 

Acceleration  caused  by  gravitation  of  the 
earth 

I 

0 

* -4^-  nondimensionalized  constant  of 

gravitation 

H 

(ft .lb) 

Kinetic  energy  of  the  gondola 

Iy>  I. 

(ft. lb. sec2 

) Gondola  moment  of  inertia  about  its  lateral 
and  vertical  axis  respectively 

1 

(ft) 

Length  of  gondola 

t 

(ft) 

Length  of  cable 
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m lb  sec^/ft  Mass  of  the  gondola 


M 

(0) 

- 1 + ^3L*z 

x 3W 

N 

(0) 

1 + 

-i  + 4 

qo> 

qi» 

q2>  q3»  q4 

(rad) 

Angles  determining  the  position  of  the  cable 
and  the  gondola  in  space  (see  Fig.  4) 

Qi» 

Q2, 

Q3 » Q4 

ft  lb 

Generalized  forces  in  the  Lagrangian 
equation  of  motion 

R 

ft 

Length  of  the  whirling  arm 

t 

sec 

Time 

t' 

0 

■ nondimensionalized  time 

T 

sec 

Turning  period  of  the  whirling  arm 

Tl> 

t2 

sec 

Period  of  the  and  q2  mode  respectively 

V 

ft/sec 

Velocity  of  the  gondola 

W 

lb 

Weight  of  the  gondola 

X 

ft 

Coordinate  in  the  direction  of  the  cable 

y 

ft 

Coordinate  perpendicular  to  the  direction 
of  the  cable. 

oC 

rad. 

Angle  of  attack  of  the  gondola 

y* 

rad. 

Angle  of  yaw  of  the  gondola 

cT 

ft 

Deviation  of  the  gondola  path  from  steady 
state  motion  measured  in  a projection  on  a 
horizontal  plane. 

0 

Damping  factor,  real  part  of  the  complex 
frequency,  positive  for  stable  motion 

TV 

0 

Length  scale  factor  of  the  model 

0 

* ma3S  fact°r 
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vk 

0 s 

- complex  frequency 

lb  sec2 
ft4 

Mass  density  of  the  air 

A 

lb  sec2 
ft2 

Cable  mass  per  length 

A 

0 

* nondimensionalized  cable  mass 

T 

0 

■ /O&V  time  factor 

CJ 

rad 

Angular  velocity  of  the  whirling  arm 

sec 

(A? 

0 

- t&'C  nondimensionalized  angular  velocity 

of  the  whirl  arm 

® 

ft2 

Cross  section  area  of  the  gaondola 

f — 

Prime  ( ) indicate  derivatives  with  respect  to  "t 

far  ( ) above  a letter  indicates  either  a nondimensional  magnitude 
see  equation  32)  or  a vector. 

Subscript  zero  (0)  indicates  a magnitude  belonging  to  the  steady  state 
condition. 
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INTRODUCTION 


When  the  Service  Branch  of  the  Equipment  Laboratory,  WADC, 
operated  the  Parachute  Test  Tower  at  El  Centro,  California  for  the 
first  time,  a serious  instability  was  observed  at  a speed  of  about 
50  mph,  that  is  about  l/lO  of  the  design  speed.  This  instability 
prevented  further  use  of  the  equipment  until  the  stability  of  the 
gondola  could  be  improved* 

Stability  difficulties  had  not  been  expected.  The  test  equip- 
ment had  been  built  on  the  assumption  that  it  would  behave  like  the 
parachute  test  tower  owned  by  Pioneer  Parachute  Co.,  Inc.  which  had 
never  shown  any  difficulties  resulting  from  instability.  The  two 
parachute  test  towers,  that  at  El  Centro  and  that  of  Pioneer  Parachute 
Co.,  however,  are  remarkably  different  in  size  and  are  not  dynamically 
similar  to  each  other. 

The  motion  performed  by  the  gondola  after  a disturbance  had  not 
yet  been  studied  either  by  systematic  tests  or  by  analytical  investi- 
gations. It  was  obvious,  however,  that  the  motion  of  the  gondola  had 
a certain  similarity  to  the  motion  of  airplanes.  Therefore  the  Air- 
craft Laboratory,  WADC,  was  requested  to  perform  a stability  investi- 
gation. 


THE  EQUIPMENT 


The  gondola  which  will  be  investigated  stabilitywise  houses  a 
parachute  and  a dummy  of  a man.  The  gondola  is  designed  to  be  whirled 
around  a vertical  axis  with  speeds  up  to  500  mph.  The  parachute  and 
the  dumny  are  released  when  the  desired  speed  is  reached. 

The  dimensions  of  this  parachute  testing  equipment  are  given  in 
Fig.  1 and  2.  The  whirling  am  has  a length  of  56  feet  and  extends 
horizontally  120  ft  above  the  ground.  The  am  is  driven  by  an  electric 
motor  of  2800  HP  controlled  by  the  Ward  Leonard  system.  Between  the 
shaft  of  the  driving  motor  and  the  axis  of  the  whirling  am  there  is  a 
gear  reduction  of  22:1.  This  high  gear  reduction  makes  the  whirl  am 
self  locking.  A disturbance  of  the  gondola  will  not  affect  the  angular 
velocity  of  the  am. 

The  gondola  has  a diameter  of  28”  and  is  112"  long  (see  Fig.  2). 

The  center  of  gravity  of  the  gondola  is  located  41^  of  its  length  behind 
its  nose.  Static  longitudinal  stability  and  weathercock  stability  is 
secured  by  tailplanes.  The  gondola  is  open  on  its  bottom  which  increases 
the  drag  of  the  gondola  considerably.  The  weight  of  the  loaded  gondola 
is  W,  - 685  lbs,  its  empty  weight  is  We  » 405  lbs* 
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NACA  0020 -G4  AIRFOIL 
CHORD  LENGTH  = 4.25' 


5CAL Ej  jo  SIZE 


The  gondola  is  suspended  from  the  end  of  the  whirl  arm  by  one 
or  more  cables  of  116  feet  length.  During  the  runs  showing  the  instability, 
a six  cable  suspension  was  used.  These  6 cables  were  housed  in  two  stream- 
lined sheathings,  each  including  3 cables.  The  profile  of  the  sheathings 
was  NACA  0020-63.  Fig.  3 shows  a cross  section  through  the  sheathing  and 
the  cables.  The  two  cable  units  were  attached  19  inches  in  front  of  and 
to  the  rear  of  the  center  of  gravity  of  the  gondola. 


THE  PROBLEM 


A solid  body  moving  freely  in  space  has  6 degrees  of  freedom.  In 
the  present  caso  of  the  gondola,  one  degree  of  freedom  is  eliminated  by 
the  cable,  if  one  assumes  that  the  cable  is  always  kept  taut.  This 
assumption  is  made  throughout  the  following  investigations* 

The  remaining  five  degrees  of  freedom  are  coupled  with  each  other 
in  a manner  remarkably  different  from  the  free  flight  conditions*  There- 
fore the  experience  gained  on  free  flying  aircrafts  and  models  cannot  be 
applied  to  the  stability  problem  of  the  gondola.  The  equations  of  motion 
must  be  adapted  to  the  present  case  and  solved  anew. 

It  is  the  purpose  of  the  investigation  at  hand  to  analyze  the  motion, 
to  find  out  which  conditions  might  cause  instability  and  which  precautions 
prevent  instability,  and  finally  to  check  the  theoretical  results 
qualitatively  by  model  tests. 

The  five  rigid  body  degrees  of  freedom  of  the  gondola  may  be 
described  by  five  angles  shown,  in  Fig.  4. 


Q ■+  - the  angle  of  the  cable  me  as  tired  in  a radial  plane, 

that  is  a vertical  plane  including  the  a raj  <3^* 
is  the  component  which  belongs  to  the  steady 
state  conditions.  It  can  be  calculated  as  a 
function  either  of  the  circumferential  velocity 
or  of  the  angular  velocity,  of  the  whirl  arm. 

^>2  • the  angle  between  the  cable  and  the  radial  plane, 

measured  in  a tangential  plane.  Precisely  speak- 
ing there  is  also  a steady  state  component  due 
to  the  drag  of  the  cable  and  the  gondola,  40 
but  <^.,0  is  small  and  will  be  neglected  in  this 
calculation* 


■ the  angle  of  pitch  of  the  gondola  measured  in  the 
tangential  plane  from  the  horizon.  If  the  gondola 
is  trimmed  accurately  the  steady  state  value 
is  zero* 
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Fig.  4a  - The  "Tangential  Plane"  and  the  "Radial  Plane" 


The  "Radial  Plane"  is  a vertical  plane  including  the  arm  OP. 

The  "Tangential  Plane"  is  a plane  containing  the  end  point  P of  the 
arm  and  the  tangent  of  the  undisturbed*  flight  path. 
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• the  angle  between  the  gondola  axis  and  the  tangent 
on  the  undisturbed  flight  path  measured  in  a 
horizontal  plane.  If  the  point  of  attachment  is 
different  from  the  center  of  gravity  of  the  gondola 
the  centrifugal  force  will  produce  a permanent 
angle  of  yaw  ^-40*  same  W be  caused  by 
unsymmetry  of  the  model  or  by  rudder  deflection. 

In  the  calculations  40  is  assumed  to  be  zero. 

• the  angle  of  roll  measured  from  the  position  where 
the  c.g.  of  the  gondola  is  in  line  with  the  direction 
of  the  cable-end- tangent. 

Precisely  speaking,  these  are  not  all  the  degrees  of  freedom  since 
the  cable  has  the  possibility  of  moving  with  an  infinite  number  of  mode 
shape 8.  These  motions  can  be  split  up  into  two  components  located  in 
the  directions  of  and  ^2* 

The  basic  mode  shape  is  very  similar  to  a straight  line  represented 
by  the  cable  with  an  amplitude  ^ 1 and  ^ 2 ♦ The  higher  mode  shapes 

have  nodes  and  will  not  be  included  in  the  present  analysis.  They  may 
be  added  later  in  a more  extended  investigation,  if  this  is  required. 

A decision  had  to  be  made  about  the  type  of  cable  suspension  to  be 
investigated.  The  tandem  arrangement  of  the  cable  units  was  suspected 
from  the  beginning  as  a possible  cause  of  the  observed  instability. 
Therefore  a mathematical  verification  of  the  instability  might  be  of 
interest.  The  following  reasons,  however,  prevented  the  investigation 
of  the  double  suspension.  The  cable  sheathing  was  so  designed  that  the 
aerodynamic  center  of  pressure  was  located  behind  the  center  of  gravity 
of  the  sheathing  and  the  included  cables.  Consequently  the  cable  sheathing 
would  be  twisted  and  high  crosswise  aerodynamic  forces  would  be  produced 
on  the  sheathing  which  could  not  be  neglected.  A sufficiently  accurate 
treatment  of  these  forces,  however,  would  not  be  possible,  since  the 
irregular  initial  distribution  of  the  twist  angle  along  the  cables  was 
not  known  and  no  data  were  available  about  the  friction  between  the 
different  assembly  units  of  sheathing.  In  addition,  the  mutual  aerodynamic 
interaction  between  the  sheathing  of  the  two  cable  units  would  have 
complicated  the  calculations.  The  results,  moreover,  would  have  been  of 
academic  interest  only,  since  it  had  already  been  decided  to  use  in  the 
future  a single  cable  suspension.  Therefore,  this  investigation  is 
confined  to  the  single  cable  suspension. 
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THEORETICAL  INVESTIGATIONS 


Performance  Considerations 


The  stability  investigation  will  show  that  the  damping  of  a 
gondola  disturbance  is  proportional  to  the  drag  coefficient  of  the 
gondola  and  cable 0 Consequently  the  highest  drag  would  provide  the  best 
stability.  Therefore  it  is  of  interest  to  find  the  highest  drag  coefficient 
compatible  with  the  required  speed  of  500  mph. 

For  lower  speeds  the  cable  angle  depends  upon  the  speed.  Conse- 
quently in  the  low  speed  range  an  exact  performance  calculation  becomes 
involved.  Fortunately,  only  high  speed  cases  are  of  interest  for  per- 
formance considerations.  For  these  cases,  however,  the  cable  can  be 
assumed  stretched  horizontally.  This  assumption  has  been  made  in  the 
following  calculations. 

The  total  drag  of  the  rotating  system  is  composed  of  the  drag  of 
the  gondola,  the  suspension  cable,  the  pipes  and  wires  composing  the 
whirl  arm,  and  its  supports. 

All  drag  components  are  assumed  in  the  fom 

3>n  * It  /Wd  cd*) 

where  /°  is  the  mass  density  of  air,  Cj)  the  local  drag  coefficient, 
which  is.  assumed  independent  from  Reynolds  and  Mach  number,  x is  the 
distance  of  a point  from  the  center  of  rotation,  d is  the  diameter  of 
the  part,  CO  is  the  angular  velocity  of  the  whirl  arm.  The  dimensions 
of  the  framework  have  been  substituted  as  given  in  the  stress  analysis 
of  the  tower.  The  gondola  suspension  has  been  assumed  as  a single  cable 
of  1-7/8  inch  diameter.  The  drag  coefficient  of  all  cables  and  steel 
tubes  has  been  assumed  to  be  1.25  based  on  their  frontal  area.  The  drag 
coefficient  of  the  gondola  is  assumed  to  be  0.2  based  on  its  main  cross 
section. 

The  result  of  these  calculations  are  plotted  in  Fig,  5.  The  upper 
curve  refers  to  the  unsheathed  cable,  the  lower  curve  to  a sheathed  cable 
with  the  drag  coefficient  Cp  - ,24  based  on  the  cable  cross  section. 

The  drag  of  the  unsheathed  cable  causes  07%  and  the  gondola  6-l/2% 
of  the  total  required  power. 
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The  lower  drag  coefficient  of 


Triangular  Sheathing 


MASS  BALANCE 


Wing  Profile  Sheathing 

Fig,  6 Forms  of  dynamically  stable 
Sheathings 


Cp  » .24  can  be  realized  by  using  a 
triangular  sheathing  as  shown  in 
Fig.  6.  This  type  of  sheatiiing  is 
considered  cheap  in  manufacturing 
cost  and  upkeep. 


A still  lower  drag  coefficient  can 
be  obtained  by  using  a regular 
symmetrical  wing  profile.  This, 
however,  will  require  a mass 
balance  and  is  considered  more 
expensive  than  the  triangular 
sheathing. 


The  available  power  is  also  shown  in  Fig.  5*  The  output  of  the 
driving  motor  is  2800  HP  and  the  maximum  power  on  the  whirling  Am  is  2240 
HP,  if  a loss  of  2(#  in  the  total  driving  mechanism  is  assumed. 

The  graph  indicates  that  the  required  gondola  speed  of  500  mph 
can  be  obtained,  if  the  drag  coefficient  of  the  cable  is  not  higher  than 

*24* 


Shape  of  the  Cable  Under  Steady  Motion 


The  knowledge  of  the  shape  of  the  cable  under  steady  state  con- 
ditions is  important.  Only  if  this  shape  is  known  can  reasonable 
approximations  be  found  and  justified.  To  begin  with,  the  shape  of  the 
cable  in  the  radial  plane  will  be  calculated.  This  is  determined  by 
the  weight  and  the  centrifugal  forces* 

In  order  to  obtain  simple  conditions,  this  problem  is  solved  in 
two  steps.  First  the  equilibrium  position  of  the  cable  is  calculated 
under  the  assumption  that  it  is  a straight  line.  Secondly  the  deviation 
from  the  straight  line  is  calculated  and  it  is  shown  that  this  deviation 
is  small* 


10 
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ion  under 


Lon  of  a 


The  coordinate  along  the  cable 
is  called  x . The  oxigin  of 
the  coordinates  is  put  at  the 

/fc(R+*3in^afdx.  «nd  P°int  p of  the  whirl  arm* 

\ Equilibrium  is  obtained,  if  the 

sum  of  the  moments  about  point  P 
vanishes.  This  condition  gives 
the  equation 

/TtffR+^sinc^oOi* 


Fig.  7 Cable  forces  for  the 

assumption  of  a straight 
cable  shape 


+ CR+%  sm  coS 

= /rwyc^Sin^+ 


■ ~ (R+gfein^o) •+  ^ + ^S\n^e]  O?  « g $ 


(R+^5inq,0)  + ^ vS/TT 


The  numerical  relation  between  CO  and  ^-0as  obtained  by  evaluation 
of  equation  (2)  has  been  plotted  in  Fig.  8.  One  sees  as  long  as  the 
angular  velocity  is  small  the  angle  Q „ changes  rapidly  with  CO.  At 
higher  CO  however,  Cj.0  changes  very  little  and  approaches  90c 
asymptotically. 
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b.  Deviation  of  cable  shape  from  straight  line 


Fig.  9 


In  order  to  designate  the  devia- 
tion of  the  cable  shape  from  the 
straight  line,  the  coordinate  y 
is  used* 

The  coordinate  x is  the  same  as 
before  (Fig.  7).  It  extends  from 
the  end  point,  P,  of  the  whirling 
arm  to  the  center  of  the  gondola. 


co5 


Coordinate  system  for 
curved  cable  shape 


Since  the  x axis  has  already  the  average  direction  of  the  cable, 
the  cable  slope  dy  will  be  small.  Consequently  the  length  element 


is 


a* 


ds  =V  d%2  + dy*  d% 


The  components  of  the  cable  tension  in  x and  y directions  are 


^9  2.  2.  _ 

. = Rx-^|  sinc^j^o  sm^o 

+ [#ri  + (£-%)/?c}  ^ cos  %0 

Ty  =[^(  R+^5mq0)+/2  (R<£+  R%-|rsin<^}axx>s^.0 


Remembering  that  {/m  ( R> ^0)j  cozCOi> 


{/m  = O 
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equation  (4)  can  also  be  written 

^ st n g.0)  co2  cos  c^e  - g sin  g:0 
-(rx  + ~ sin  g.e)  co2- cos  + xg  smcj. 


Tv  * /°c 


In  this  form  it  is  ^asily  seen  that  Ty  changes  its  sign  in  the 
neighborhood  of  , which  indicates  that  the  cable  curve  can  not 

deviate  too  much  from  the  straight  line. 

The  differential  equation  of  the  cable  curve  is 

dt 


= J*L 
dX  Ty 


(5) 


COS1 

which  gives  after  factoring  "3771 
* 


nj  « Co-t  3. 


fri( +/?c  ^ 51  n 5 £2/7 

-/^c(fR^-gt&A?9e)x-  ^602-5in^-0  x2- 

C R+ JC.  s\  n C R <.1-  1 (n-n  3 Co-6  ^.o  76 

-/l(R  co*-+  g co-t  4Jx-  ^ ca>\si  n cj.b  x/ 


The  term  in  the  denominator  is  removed  by  dividing  by  the  numerator: 

g(ixui$0+C  crfcfe)  , 

+ /g(R*. + n^o)}coz  ^ 

+ cot$6  -/£[Ra>Vg  cot^0) X 

- ^ COZStn  ^-0X2  _ 

After  splitting  the  fraction  up  into  linear  partial  fractions  and 
integrating,  one  obtains 

tf  = co tg^t  cotq.^  jj^  W,_g-)  + K£lnO“^)]} 

(6) 
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ct+<*  b 

oL~/S 


K1  - 


*2 


CL  +/g  b 
oL~/g 


a - — {/M  +/%sC) 


A - [/w(Rtiesinq.&)+/g(f?5e+  |^%m^}6o2--+('/W+/^^;^  cotCj.0 

b » ^<3C 


B « /?c  [RcoN-g 
m OO2  sin  g.. 


The  shape  of  the  cable  is  plotted  in  Fig.  lOlbr  q0  ■ 45*, 
the  angle  at  which  the  maximum  curvature  of  the  cable  can  be 
expected.  One  sees  that  the  deviation  from  the  straight  line  is 
less  than  1%  of  its  length.  The  substitution  of  the  cable  shape 
by  a straight  line  is  thus  justified. 
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The  Pendulum  Modes  Treated  As  Uncoupled  Degrees  of  Freedom 


a.  Outward-Inward  Oscillations : This  motion  is  restricted  to  the 

radial  plane  (Fig.  4a).  The  forces  acting  are  the  weight,  the  centrifugal 
forces  and  the  cross  forces  on  the  gondola*  Assuming  that  the  gondola  has 
a high  enough  weathercock  stability,  the  gondola  will  turn  its  nose  always 
into  the  resultant  wind  so  that  the  cross  forces  become  zero*  Therefore 
the  airforces  will  be  neglected  in  the  treatment  of  the  uncoupled  motion* 

The  total  moment  about  the  point  of  suspension,  P,  is 

M * £rt»[R+£sin  (q.c+  q.,)]  £ + 3- ^cos 

- & %-} 4-sin (q.0+q.,) 


The  equation  of  motion 


£ 


q.  - g,  - o 


becomes  4 

+/^(Ry  + ^5in %V}  6oz5 1 n qc-  ^ } co5*q.0 0zjq.*0 


or 


9-, 


H- 


'/>n  +/zf. 


COS 


The  related  frequency  is 


**  V2fr 


® OOf 


vm 


11 1 s'^->  US*  COs\ 


-'wl  +/%lir 


J 


CO Sg. 
3C  602- 


(8) 


The  period  Tltl  « -^r-  is  plotted  versus  the  whirling  frequency  in 
Fie.  11,  where  it  is  marked  Tia  . 

The  period  of  the  uncoupled  outward-inward  pendulum  motion  reaches 
about  90$  of  the  whirling  period. 

b.  Fore  and  aft  oscillations:  This  type  of  oscillation  occurs 

in  the  tangential  plane.  The  restoring  forces  are  again  given  by  the 
components  of  forces  related  to  the  gravitational  and  the  centrifugal  forces. 
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Coriolis  forces  do  not  enter  the  investigation  of  the  uncoupled 
oscillations,  since  these  are  typical  coupling  forces. 


The  restoring  moment  is  produced  by  the  gravitational  forces  and 
by  the  centrifugal  forces.  The  centrifugal  forces  act  horizontally  and 
point  in  a direction  away  from  the  axis  of  rotation.  The  centrifugal 
forces  can  be  resolved  into  two  components,  one  producing  a moment  tending 
to  restore  0-2  > to  zero,  the  other  producing  a cable  force  and  a moment 
which  trie«  to  increase  C^q.  (See  Fig*  13) 


The  restoring  moment  produced  by  the  centrifugal  force  on  the 
gondola  is:  cm  £r+  a/'Stn^eC  where  /&'  is  the  angle 

defined  later  in  Fig.  14  and  in  equation  15.  Substituting  the  value 
for  sin  and  adding  the  contribution  of  the  cable  one  obtains  the 
total  restoring  moment:  -p 

CQn 

Sin  <3.0 

"o  o 


=-[- 


+ C05q.o  d-il  q . 

' "X  -1 


= - )( 


sin%, 


-t-  $cosq.a£.) 


The  equation  of  motion  is 


(ar\X.* 


3 


) ^ +■  (> 


2. 


which  has  the  frequency 


(9) 


The  period  TSu_  =•  is  also  plotted  in  Fig.  11.  For^.4=5j  *1*2*. 

starts  out  with  less  than  1/3  of  the  value  of  T and  becomes  higher 
than  T for  greater  than  60°. 

The  fact  that  all  three  periods  Tq  , T2,  and  T are  of  the  same  order 
of  magnitude  leads  one  to  expect  strong  couplings.  This  will  make  the 
evaluation  of  tests  difficult  unless  one  of  the  two  modes  is  much  more 
damped  than  the  other.  It  also  makes  an  investigation  of  the  coupled 
motion  necessary. 
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Derivation  of  the  Equations  of  Coupled  Motion 


The  investigation  of  the  stability  of  the  steady  state  motions 
of  the  gondola  is  the  main  purpose  of  the  present  report.  This  will 
be  done  on  the  basis  of  the  theory  of  vibrations  developed  by  Lagrange. 

The  degrees  of  freedom  considered  are  described  by  the  coordinates 
^1,0-  which  have  been  introduced  before.  This  means  that  the 

rolling  motions  of  the  gondola  and  the  torsional  and  translational 
motion  of  the  cable  in  itself,  which  do  not  affect  the  position  of  the 
cable  ends,  are  neglected. 


In  order  to  eliminate  the  internal  forces  produced  by  the  cable, 
the  equations  of  motion  will  be  used  in  the  Lagrangian  form.  Using  the 
symbol  H for  the  kinetic  energy  and  Q^,  Q2>  Q4  for  the  generalized 
forces  belonging  to  the  degrees  of  freedom^,  ^-3 * ^4*  respectively, 
the  Lagrangian  equation  of  motion  can  be  written: 

cl  an.  _ sh  _ 

dt  a q_,  $<q , 

d ££L  _ 3H  _ n 

dt  o)  2-  . 

d 3H  3 H _ 0 (10) 

dt  aq3  ^3 


d_  2H.  _ aid 

d't  Q 3^. 


0-4- 


The  kinetic-energy  will  be  calculated  in  the  rotating  system. 
Consequently  the  centrifugal  force  and  the  Coriolis  force  will  be  included 
in  the  generalized  forces  Qn. 


The  kinetic  energy  H can  be  split  up  into  two  parts,  the  kinetic 
energy  Hg  of  the  gondola  and  the  kinetic  energy  Hc  of  the  cable. 

The  velocity  has  two  components  each  one  in  the  direction  of  the 
coordinates and  C^2*  For  the  gondola  these  components  are 

Vig  = -i  4, 

v&g  - 
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For  a mass  element  of  the  cable  having  a distance  x from  point 
P the  equations  are: 

V,c  = ^4. 

V2C  ~ 

Consequently 


H - H3*Hc  = — 


2 


+ =:  q: 


+ ¥ ^4 


The  generalized  force  Qn  is  the  work  done  by  the  external  forces 
during  a virtual  displacement  divided  by  the  displacement  . 

The  external  forces  are  the  gravitational  forces,  the  air  forces,  the 
centrifugal  forces,  and  the  Coriolis  forces,  that  is  all  forces  except 
the  cable  force,  which  is  an  internal  force* 

The  direction  of  the  gravitational  force  is  vertically  downward. 

The  centrifugal  force  is  in  a horizontal  direction  outward.  The  Coriolis 
force  can  be  written  in  vector  symbols 


(12) 


and  may  be  split  up  into  components  according  to  the  two  components  of 
the  velocity  V . Since  all  vectors  of  the  form  £5*  <3-  are  perpendicular 
to  the  vector  <-o  which  is  directed  vertically  downward, ^ they  are  all 
contained  in  a horizontal  plane.  The  component  acoss^  c0-R-( points  back- 
wards, while  the  component  £ JCtoij,  points  away  from  the  axis  of  the  tower, 
(See  Fig.  12)  * 

The  lift  force  L on  the  gondola  rolls  with  the  gondola  and  has 
almost  the  direction  of  the  cable  (See  Fig.  12).  Drag  and  cross  forces 
are  perpendicular  to  the  lift. 

It  is  assumed  that  the  gondola  is  trimmed  for  Cl  * 0 and  Cc  = 0 
and  that  it  moves  at  minimum  drag  coefficient  as  long  as  the  speed  is 
steady.  This  means  that  ^^=0  and  O • 
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I 


Centrifugal  fo  rce : — /??7  CO  x ( 60  x r) 
Coriolis  force:  -2 ./ni  oo  KY 
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Subscript  p means  the 
magnitude  is  not  shown  in 
full  value,  but  in  pro- 
jection only. 


Consequently 

Gt,  =^^2[R+£5im(cfo+^]<cos(^0+^i)  -/wgr^sm^+q.) 
rjf* 

+ co2- J Cos(^<^) <J*-/g 9 S\n  (<^q.,) J ** d*. 

+ (CD  stn^-Q.)  ^ V£®^  + £/wi^2go  GosC^.o+9-,)^^ 

+ e/Z^c'osC^+q.,) y ; ^zJ^q.z 


or 


Q.,  * /w{  £0*  [R+<sin(^^)J^cos(^o+^  -g^sm^+^J 

+ ^ 6o2  [ R 1%  J3s \n  C^0+^)] cos (c^q. ,)  2sm  (^B+%)f 

+ (Cd-^C)/5^  V2®^.  -+  COS^-t-^W^dJ) 


Before  the  generalized  force  Q2  is  derived,  the  location  of  the  - 
centrifugal  forcefon  'toe  gondola  relative  to  the  virtual  displacement  0^2. 
will  be  described,  tfe  is  located  in  a horizontal  plane  and  points  in 
the  direction  from  the  axis  of  rotation  outward  (See  Fig.  12).  It  is 
resolved  into  three  components,  ^-^has  the  direction  of  the  cable,  ^iis 
perpendicular  to  the  cable  in  the  direction  of  increasing  Q_.  , and 
is  horizontal  pointing  in  the  direction  of  decreasing  • 

One  reads  from  the  rectangular  spherical  triangle  Fig.  13. 

cos  p'g  = c.osq.2.  cos  (q.a+  qp 

.3 in  Ql  = 

t-2  sm  p'G»  - coszq.zco£2^q.0+^.iy 


When  q.c  approaches  zero  and  if  and  are  small  one  obtains  for 
^2  ^e  limit  value  ^/z  a3  one  would  expect. 

For  finite  values  of  ^-D  and  small  values  of  <9..  and  q.2  on®  can 
assume  x ' 


si  n %o 


(14) 
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Hence,  when  is  neglected  in  comparison  to  Q_o  : 


. R Q-a 

LR+sCsm^.-i-q.,)]  smqB 


(15) 


Using  this  result,  the  component 

becomes 

^3 l3  =•  *6$^* 


of  the  centrifugal  force 

_ sm  R c.o'2,  „ 

- St 


This  force  is  perpendicular  to  PG".  Similarly  one  obtains  for  the 
corresponding  component  of  the  centrifugal  force  on  a cable  element 


4 


Z- 


After  adding  the  contribution  of  gravitation,  of  the  Coriolis  force, 
and  of  the  air  forces  one  obtains  the  generalized  force 


£ 


r i—  _pz-_  • 1 

Q.z  - > ■ gw  £ 

- 2.  9- cos  C /w+y#  - x>«£ 


(16) 


This  equation  is  not  valid  for  = 0,  since  equation  (14)  has  been 
used.  The  case  9-e*  however,  is  of  no  interest.  The  contribution  of 
the  lift  force  is  but  small  and  of  second  order  and  has  been  neglected. 

The  Coriolis  force  might  be  treated  in  the  same  way  as  the  centrifugal 
force.  Since,  however,  the  Coriolis  force  i3  small  but  of  the  first  order, 
the  moment  which  originates  from  it  by  the  multiplication  by  9- a would  be 
small  and  of  the  second  order.  This  component  is  therefore  neglected.  The 
generalized  forces  and  include  aerodynamic  moments  only. 

+ a?) 

<34  - 4.)  ® 2 -f  V2  (18) 
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The  mass  forces  do  not  contribute  to  the  moments  Q3  and  since 
Qo  and  Q,  are  rotations  about  the  center  of  gravity.  When  the  values 
(17)  and4(18)  are  substituted  into  the  equations  of  motion  it  will  be 
considered  that 


(19) 

“ ^-4  s 7 R + arSnq.0-)6o 

(20) 

These  values  are  now  substituted  into  the  Lagrangian  equations  of 
motion  (10). 


d. 

d-t  ^ 

H 


3C3-  " 


ri 


= O 


ct_ 

d-t  <3c|3 
3 ^3 


ft  **%»+$*) 

sm 


= o 


c^e: 

£>H 


g)H 


^4 


Consequently  the  equations  of  motion  are: 

('W+/9.  = 


- 

= O 

- 1*  4* 

= o 


'Tn [\r +£ sm(q.0+  % )\  ^ cos (^Q+  co*- g^s\ n (<$0 ,+q.,') v 2. £2to  cos  <f<q0+^() 

•1-/?  r[R^|3s\n(q5v%^a)^s(v%l')  - J o)co5(^-q.l>)q^ 
+ ®'^C<3-4',r  CR+^sm^.6)oo) 


(21) 
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(22) 


(*»  + + 2 (*»  4-/2  j)  4-,  CoS  % «> 

■+[gcos^0+q.,)  JW+/2  f*]  R*  <^}  3-z 

+ /~VZ'CJ>  ® £ « O 


(23) 
O (24) 


®2|V^  <=  O 


(R+5fSm^)6y)  + lq^  3-4-1  ® 2 2^  V 


^4 


Linearization  and  Nondimensionalization 


of  the  Equations  of  Motion 


The  investigation  will  be  confined  to  small  disturbances  to  the 
steady  motion.  Therefore  it  is  permissible  to  neglect  terms  which  are 
proportional  to  quadratic  and  higher  powers  of  the  disturbances.  If  the 
angle,  by  which  the  gondola  lags  the  whirling  arm,  is  disregarded,  the 
disturbances  are  On,  2>  and  Cj_4.  Terms,  which  do  not  include 

these  variables  belong  to  the  steady  state  motion,  which  has  already  been 
taken  care  of.  Consequently  only  those  terms  are  retained  in  the  equations 
of  motion  which  are  linear  inC^i,  Q.  2,  <^.3,  and  <3.4. 

If  one  uses  the  abbreviation: 

VQ  - (r  + sCsi tog..)  CO 
the  square  of  the  speed  can  be  linearized  as  follows : 


V2-  a \r0z+  [R+^svnq.0][cosq.oafg|+  60  q.2]  (25) 

After  linearization  one  obtains  from  equations  (21)  to  (24) 


[Ri-^s\nq.^afswq.0^,+5C2’cos^06o^.l  - g£coscf0q.,  -t-2.£zcos  ^eu;g3 

+/*[-[■ 


^Cc 


+ g-'cos^^j. , - grfcos  , + zfcos  qbco  <jj 


?I26) 


+ (C;D  3/3  ) (R  + -jLti\r)  cp  uj]  2 
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* *.  . C*L.  \ 

/Jr,  “ ■ 


’%_+  aT^OCf,  C.05^, 

01 l±&i. 


--  i 3 CP~  O-o  (,w£+/i  ) 

^ <C 

^/^(R  + ^Gin^XcQS^CO^,^  «^)  Cp<X>^. 


O 


(27) 


^.-C  §&“%  ■’■ 


aC 


2Z^ 


Q 

T3 


) 


® i 


/>  2. 
p Vo 


= o 


(28) 


(29) 


These  equations  will  be  transformed  into  the  nondimensional  form. 
The  nondimensional  form  has  the  advantage  of  a direct  comparison  between 
full  scale  and  model  tests  and  making  clear  which  combinations  of  the 
parameters  are  essential.  Moreover,  the  coefficients  of  the  equations 
become  all  of  the  same  order  of  magnitude  which  makes  it  more  convenient 
to  handle. 

For  this  purpose,  the  following  substitutions  are  made:  Instead 

of  seconds,  the  value  7r  Is  used  as  a unit  of  time,  which  is  defined 
by: 


The  mass  factor: 


77 


^77-7 

/o  CS>  V. 


/m 

/3<&l 


and  the  following  abbreviations  are  introduced 

-*5  -J-  z' 

£L  7-2- 

Jbz  - JL  JXi 

/g  cC  ^ ' Z. 

/m 

4 M - H'r  - 


<70  = 

g = 

/c  = 

iT  - 
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3/m 


(30) 


(31) 


(32) 


Primes  ( f ) are  used  for  designating  derivatives  with  respect  to  t. 

After  multiplying  equations  (26)  and  (27)  with  and  after 

dividing  equations  (28)  and  (29)  with  ^Vo1®  l one  obtains 

M off  + ( g-+5in^0)£o2^m  ^ - cos^or^  ,+•  § cosq^-zcosyo^.' 

1 cos  ^4r|co^0“j 
* z 'k^ (C;p “ 1^0  ( ^ ) = ° (33) 


Mq[^  + 2MO)  g.,  cos  q.B 

_ 

4 


[5  <*» *.  C I t 4^  £ “*}  ^ 

+Jr(  J + 7 •sin^.)(cosq.e60^,+  604',I)  Cj, 


= o 


(34) 


^ -H 


\ / ^Cm. 

of  ®C 


% 


= o (35) 


J_ 

£ 


^S-I  \ _v  Q ' 

{'R+^Stng,,)  CO  J r0C^  ^"4 


- o 


(36) 


After  collecting  terms  the  following  equations  are  obtained: 

M < - **  (c»  - Wt)  (R 

~ +.Sing.e4-  ^ -v^-  s\n<^)j  co' zst to  ^ 6 - aJ^o^, 

+/^(f  55*COS-%^)  CO5^6  + §CoS0^A  Qh 

- 2 M COS<^.b  OOc£_i  |^CJ>  - l^f)  S-4  = O 


ac.osq.D  M co 

4 hAC^  4 


9.; 

Jc  (f+-?sm^*i  <3JC»44 

+ {(3cos^.H-3Ss_ro*Xl-r't)}^ 


=.  o 


(38) 


(39) 
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Jk 


2 

a 9-3 


1 3 C >>v 

2 


-L 

2 3 ^ 


9-3 


= o 


(40) 


„!  iCn  £ ^q/.  JLV  ..laCn 

2 9/5  ('R^smq^aj)  n ^^4  2 9^.4 


1 

2 9/5 


(41) 


The  solution  of  the  equations  is  assumed  in  the  form 


0^  * n — An  0. 


■£-$,+ *>*>4 


(42) 


If  the  solution  is  not  identically  zero,  which  would  represent  the 
case  of  undisturbed  motion,  which  is  of  no  interest  here,  the  determinant 
of  the  system  of  equations  (38)  to  (41)  must  vanish.  This  determinant 
has  the  form 

bl2_V 


(45) 


cuD2,+  b,,n)  + a, 


b^  V t 0-2\  d"  bgjf'l)  4'  Q 


Z2. 


b41i) 


<^+M+a23 


where 


Cjj  = O i ’?  ~ N4 


zz 


bi|  = + q„Vo  2 ^ _ 

= -t-  MsmqJ  a3^u',q(.-  CClM  cos*4B+  s6+^coS^- 

bl2.~  - 2 M 60  Cosq.e 

b2)  = 2M  53  cos^e  ~ --  b . 


(44) 


”^V  V.  z 


;ih.q.t)cr.-o(  co^C^  « 60  CD  coc  q_, 


Q44 

= 0 

■0^4. 
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It  follows  from  the  determinant  (43)  that  the  pendulum  degrees  of 
freedom  9-  x an<*  9-2  are  independent  from  the  rotary  degrees  of  freedom 
9-3  and  if  either  a^  or  is  zero.  can  become  exactly 

zero,  if  Cp  ■ , although  even  a slight  deviation  from  the  exact 

zero  gives  ai4  an  appreciable  value  because  of  the  large  numerical  value 
of with  which  (Cp—  ) is  multiplied. 

The  other  condition  under  which  the  separation  of  the  two  groups 
would  take  place  is  that 

H =r  — — 5>  =£• 

41  Z 

would  be  negligibly  small.  Since  is  the  order  of 

.02  and  since  C^OOis  very  likely  to  have  a Ji  type  characteristic  with 
an  almost  horizontal  tangent  for  * 0,  would  also  be  very  small. 

Therefore  the  assumption  that  k/ti  • 0 is  a gaod  approximation. 
Nevertheless  it  should  be  kept  in  mind,  that  from  here  on  the  investigation 
is  somewhat  simplified. 

If  b^  * 6,  the  determinant  (43)  splits  up  into 

C„-iA  b,,nP+GL„  bl2_-i> 

(C33t A b33i)  + a33)( C44-A  bA4'0+CLH)  “°  (45 ) 

to2  + cl  2 1 C$)  + b2£V+  C/i-i 
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Instead  of  an  eigth  degree  polynominal , two  quadratic  and  one 
quartic  equation  has  to  be  investigated. 


Stability  of  the  Coupled  Pendulum  Degrees  of  Freedom 

The  quartic  determining  the  stability  of  the  pendulum  degrees  of 
freedom  and  ^ 2 can  be  written: 

A4-a)4  + •+  A-ia)*  + A,-0+  A*  =*  O (46) 


Here,  the  following  abbreviations  are  used: 

A4  ■ Cll  c22 

a3  " C11  b22  + c22  bll 

A2  " C11  a22  + bll  b22  + C22  all  ” b21  b12 
A1  " bll  a22  + b22  all  “ b12  a21 

Ao  * all  a22 

The  conditions  for  stability  are 

A^  > 0;  > 0;  A 2 >■  0;  A^^  Oj  Aq^  0 


(47) 


'R  - A,A7  a,  - Af  A4  - Ao  A; 


O 


If  the  Routhian  discriminant,  R,  is  going  to  be  chosen  as  a measure 
of  stability,  the  coefficient  A^  of  the  frequency  polynomial  must  be  made 
equal  to  one.  This  means  that  each  coefficient  A is  divided  by  A^  * M^. 
If  we  retain  the  coefficients  as  being  given  by  equations  (47),  the 
discriminant  Routhian  has  to  be  divided  by  M°.  One  has 

M*  = [a,(AzA3  - A,  A4)  - Ao  A|  j 

Substituting  the  values  from  equations  (47)  and  (44)  and  considering 
that  b~n  “ 0 one  obtains : 

A2A3-A,  A4  - C3>MtCaM+  4M  a5^cosz^.c')-MzCc3)a,,  + aM  afcos^Cj,) 

= CjjM2, 2MwzCCS^] 


A|(AZA3— AA4)  - c£ Mt(an+  2 M coz-cos%')(a22-t  ^ m of  cos?-c^o) 
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A0Aj  - cl„  a«.  M2  Cjf 

= ^4.^C<X„4  Ct£Z)Z  M 602-  cos%  + 4 M .0-  cos4q.  o} 

-j^fc  = a3zcos%0^a,,  + cizz-v  2 M a 52  cosz^DJ 

The  values  of  aj^and  a22  are  also  substituted  according  to  equations 
(44).  The  abbreviation 


N 


(43) 


is  used.  It  is  • _ 

a„  » m{^N  co*3mSk+u>  Omz90-co5%>)+  <3  N cos^^j 

<*«.-  M M £§  cos  <*..  + 

Hence  __ 

yfi.  - co4‘costiq.^  l + M ^ os  ^(sm?l-i+sin  S°^ij 

After  substituting  the  nondimensional  values  3 and  co  by  the  original 
values  one  obtains 


R 

Mb 


■sC-t zCQ5x<su /7n^ 

/°*  (^YR-v^sinq.^^ 


£cf 


cos^0+  ^(5inq..-v 


(49) 


Since  O 4 ^ 4 ^ it  is  obvious  that  — g-  can  never  become  negative. 

Consequently  the  only  possibility  of  instability  is  that  one  of 

the  coefficients  A-^may  become  negative.  Numerical  calculations  vrill 
show,  however,  that  this  is  not  the  case.  Consequently  no  dynamic 
instabilities  are  to  be  expected. 


The  result  of  a numerical  investigation  of  the  stability  of  the 
system  with  a single  torsion  free  cable  is  reproduced  in  Table  III.  The 
calculations  are  based  on  the  data  of  Table  I.  The  calculations  are  made 
for  the  fully  loaded  gondola  of  W * 635  lbs  and  for  a cable  with  a weight 
of  5.17  lb/ft  including  the  weight  of  the  sheathing. 
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The  Table  III  lists  the  coefficients  A0,  A]_,  and  A2  of  the  frequency 
polynomial  for  the  range  of  from  5°  to  85°.  A3  and  A4  have  the 

constant  values  A3  « .2326;  A4  - 1.670.  The  drag  coefficient  of  the 
gondola  including  part  of  the  cable  drag  is  assumed  to  be  Cp  ■ 0.24. 

This  value  might  actually  be  somewhat  higher.  This  difference,  however, 
would  affect  the  stability  only  very  slightly. 

The  frequency  polynomial  has  been  solved  by  means  of  the  method 
developed  by  H.  C.  Gebelein*  which  is  considered  especially  suited  for 
cases  like  the  present  one,  where  the  roots  are  located  relatively  close 
to  each  other.  The  damping  ratios  ^ 1 and  ^ 2 are  extremely  low,  i.e,, 
in  the  order  of  about  10~7  and  10“3.  Because  of  equation  (49),  however, 
one  is  sure  that  an  inaccurate  assumption  or  a variation  of  the  data  will 
not  lead  to  instability. 

The  periods  T^  and  T2,  which  were  calculated  previously  neglecting 
the  mutual  couplings,  are  also  listed  in  Table  III.  A comparison  with  the 
former  values  shows  that  the  coupling  decreases  the  periods  at  all  cable 
angles  CJ_0  . The  cable  angle,  however,  for  which  the  two  periods  are  equal 
to  each  other,  remains  the  same  (52.4°). 


MODEL  TESTS 


The  theoretical  investigations  shown  on  the  preceding  pages  had  to 
be  made  without  reference  to  other  reports  since,  no  analytical  treatment 
of  the  same  or  a similar  problem  was  known.  Therefore  it  was  desirable 
to  verify  the  results  by  tests  at  least  qualitatively. 

Checks  on  the  full  scale  equipment  were  out  of  the  question  because 
of  required  time  and  cost.  Therefore  recourse  was  made  to  tests  on  a 
1:20  scale  model* 

It  is  obvious  that  a dynamically  similar  model  in  an  exact  sense 
could  not  be  built.  While,  for  instance,  the  full  scale  Reynolds  number 
of  the  round  unsheathed  cable  is  above  the  critical  number,  the  corresponding 
Reynolds  number  of  the  model  was  below  the  critical  number* 

These  insufficiencies,  however,  were  not  considered  serious,  since 
the  theory  had  revealed  their  secondary  importance.  The  intent  of  the 
model  tests  was  not  so  much  to  reproduce  the  characteristics  of  the  full 
scale  equipment  as  to  check  qualitative  theoretical  results  similar  to 
the  following.  The  frequency  of  the  two  pendulum  modes  is  about  equal 
and  as  long  as  the  cable  angle  is  below  70°  both  periods  are  shorter  than 


* H.  C.  Gebelein:  Treatment  of  Fourth  Degree  Equations  - AF  Tech  Rpt  No.  5992. 
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the  period  of  turning  of  the  arm.  All  configurations  going  to  be  tested 
will  be  stable,  but  the  damping  will  be  very  low. 

It  was  especially  the  latter  result  which  would  represent  a sensitive 
check.  The  only  instrumentation  needed  would  be  some  means  of  recording 
the  motion  following  a disturbance. 


Dynamic  Similarity 


The  dynamic  similarity  of  the  model  was  based  on  the  length  scale 
factor  -A.  • 1:20  and  on  the  requirement  that  the  direction  of  the 
resultant  of  centrifugal  and  gravitational  forces  should  be  the  same  for 
the  model  and  for  the  full  scale  equipment.  Since  the  gravitational 
acceleration  g is  the  same  for  model  and  full  scale  equipment,  the 
centrifugal  acceleration  had  also  to  be  the  same.  This  means  that 

the  circumferential  speeds  had  to  have  the  ratio  VUjeL  * = 1 ; V20 


.3 

Masses  were  reduced  by  the  factor  -A.  « 1:8000.  This  applies  to 
the  mass  of  the  gondola  and  of  the  cable.  The  similarity  of  the  air  forces 
is  fulfilled  automatically. 


For  two  reasons,  no  effort  was  made  to  obtain  similarity  of  the 
moments  of  inertia  of  the  gondola.  First,  the  moment  of  the  inertia  of 
the  present  gondola  was  not  known  and  the  data  on  the  new  gondola  which 
is  going  to  be  designed  were  not  yet  available.  Secondly,  the  theoretical 
investigations  had  shown  that  there  would  be  very  little  coupling  between 
the  rotational  degrees  of  freedom  and  the  pendulum  degrees  of  freedom 
which  are  of  interest  here. 


Model  Test  Equipment 

It  was  out  of  the  question  to  build  more  than  one  gondola  model. 
Since  the  full  scale  gondola  will  later  on  be  suspended  with  a single 
cable  attached  above  the  center  of  gravity,  the  single  cable  suspension 
was  reproduced  by  the  model.  The  dimensions  of  the  new  gondola  were  not 
yet  available.  . Therefore  the  model  was  built  in  the  form  of  the  old 
gondola  adapted  to  a single  point  suspension  as  being  used  at  El  Centro 
by  that  time.  The  cable  was  attached  to  a yoke,  which  connects  the  two 
former  point  of  attachments  (See  Fig.  15  and  16). 

This  arrangement  has  a large  lateral  area.  Therefore  the  condition 
CD  - ■ 0 was  very  likely  not  fulfilled.  Hence  it  was  expected  that 

tests  with  this  model  would  also  indicate  to  what  degree  a violation  of 
the  assumption  Cp  - =■  0 would  modify  the  theoretical  results. 
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YOKE  CONVERTING  TWO  POINT 
SUSPENSION  INTO  SINGLE 
POINT  SUSPENSION. 
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Fig.  15  Gondola  With  Heaviest  Weight,  Wj,  and 
Medium  Drag  Plate,  T>2 


Fig.  16  Gondola  With  Heaviest 


The  driving  mechanism  consisted,  during  the  first  tests,  of  an 
electric  motor  designed  for  high  gondola  speeds.  At  the  test  speeds, 
which  were  much  lower,  the  motion  was  erratic.  Even  when  a gear  reduction 
of  11:1  was  used  the  motion  was  not  even  enough  for  reliable  measurement. 

Therefore,  the  electric  motor  was  replaced  by  a compressed  air  motor. 
Though  this  gave  a much  smoother  drive  serious  fluctuations  of  the  speed 
were  encountered  when  other  tools  connected  with  the  same  compressed  air 
tank  were  in  operation.  The  available  pressure  regulator  could  not  cope 
with  these  pressure  fluctuations.  Finally  an  electric  motor  controlled 
by  the  Ward  Leonard  system  was  installed,  which  gave  satisfactory  smoothness 
of  the  turning  speed.  (Fig.  17) 

The  disturbance  of  the  motion  was  produced  by  a jet  of  compressed 
air.  The  jet  was  formed  by  a nozzle  mounted  on  a tripod  (Fig.  18).  A 
suitable  direction  of  the  nozzle  was  found  by  trial  before  a measurement 
was  made.  The  distance  of  the  nozzle  from  the  gondola  path  wa3  so  chosen 
that  the  disturbance  was  of  a reasonable  magnitude. 

The  purpose  of  the  tests  was  to  obtain  numerical  results  on  the 
periods  and  the  dampings  of  the  disturbed  gondola  motion.  The  easiest 
way  was  considered  to  be  the  taking  of  time  exposure  pictures  of  the 
motion.  For  this  purpose,  the  nose  of  the  gondola  was  made  of  transparent 
plastic  and  an  electric  light  bulb  was  installed  underneath  of  the  plastic 
cover  (See  Fig.  15  and  16).  Light  bulbs  were  also  mounted  on  the  end  of 
the  turning  arm  and  at  the  turning  center  (Fig.  18).  The  traces  of  these 
two  latter  light  bulbs  helped  in  aligning  and  comparing  photographs  of 
different  tests. 

Two  cameras  of  4”  x 5"  size  were  used  for  recording  the  motion.  One 
camera  was  located  above  the  model  tower.  By  means  of  a plumb  line  the 
optical  axis  of  the  camera  was  lined  up  with  the  axis  of  rotation.  The 
view  taken  by  this  camera  is  shown  in  Fig.  20.  The  other  camena  was  on  a 
tripod  facing  the  model  tower.  The  distance  of  its  lens  from  the  axis  of 
rotation  was  29  feet,  2-1/2  inches.  Its  optical  axis  had  the  same  height 
above  the  floor  as  the  turning  arm.  A view  taken  by  thi3  camera  is  given 
by  Fig.  21. 

Eighteen  different  test  combinations  of  gondola  weight,  cable  weight 
and  additional  drag  ware  tested,  as  listed  in  the  following  table. 


Run  No. 

Gondola  Weight 

Cable  No. 

Gondola  Drag 

1 

Wi 

Ci 

Dl 

2 

w2 

C1 

Cl 

3 

W3 

C1 

Dl 
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Fig.  17  The  Driving  Mechanism 

WADC  TR  52-72 


39 


WADC  TR  52-7S 


1+0 


WADC  TR  52-78 


Fig.  19  Model  Set-Up  Shewing  Use 

Of  Protractor  For  Speed  Adjustment 
l4l 


Fig.  20  View  of  Overhead  Camera 


Fig.  21  View  of  Ground  Camera  Showing 
Model  in  Operation 

1+2 


-Trace  of 
Tip  Bulb 


-Trace  of 
Gondola 


-Model  Tower 


WADC  TR  52-7S 


Run  No.  Gondola  Weight  Cable  No.  Gondola  Drag 


4 

W1 

c2 

D1 

5 

w2 

c2 

D1 

6 

W3 

c2 

D1 

7 

W1 

C1 

d2 

8 

w2 

C1 

d2 

9 

w3 

Cl 

D2 

10 

*1 

C2 

D2 

11 

«2 

C2 

d2 

12 

w3 

c2 

D; 

13 

Wl 

Cl 

14 

w2 

Cl 

d3 

15 

W3 

Cl 

d3 

16 

Wi 

C2 

d3 

17 

w2 

C2 

d3 

18 

w3 

C2 

d3 

The  numerical  values  of  the  gondola  weights,  cable  weights  and  dimensions 
of  drag  plates  are  listed  in  Table  II.  Each  run  consisted  of  measure- 
ments at  5 different  speeds  corresponding  to  cable  angles  Q0  * 30°,  40°, 
50°,  60°,  80°.  Each  speed  was  run  at  least  three  times.  Consequently 
the  total  sum  of  test  points  was  over  270. 

The  procedure  on  each  test  point  was  as  follows.  The  speed  was  so 
adjusted,  that  approximately  the  intended  cable  angle  was  obtained.  For 
this  purpose  a protractor  was  used  as  shown  in  Fig.  19.  The  accurate 
speed  was  determined  by  means  of  a stop  watch.  The  speed  measurement 
was  repeated  before  and  after  each  measurement. 

Next  the  air  nozzle  was  adjusted  30  that  a disturbance  of  suitable 
amplitude  was  obtained.  The  lights  on  the  model  were  switched  on  and 
the  lights  of  the  test  room  switched  off. 
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Turns  were  counted  from  one  to  seven.  Shortly  before  "three”  was 
pronounced,  the  airjet  was  turned  on  and  the  shutter  of  the  cameras  were 
opened.  The  air  jet  was  closed  as  soon  as  the  gondola  had  passed  it. 

The  shutters  were  left  open  up  to  "five"  or  "seven"  respectively.  The 
path  of  the  undisturbed  motion  was  recorded  on  a separate  film.  Examples 
of  the  pictures  taken  by  the  overhead  and  the  ground  camera  are  shown  in 
Fig.  22  through  27. 


Evaluation  of  the  Tests 


In  order  to  evaluate  the  records  taken  by  the  overhead  camera, 
enlarged  prints  of  the  size  14  x 14  inch  were  made  and  the  deviations  £ 
from  the  steady  state  path  were  measured  and  plotted  versus  the  azimuth 
angle.  In  cases,  where  the  shutter  was  left  open  until  "seven"  almost 
unsurmountable  difficulties  were  met  in  following  the  trace  through 
the  many  recorded  turns.  A plotting  of  7 tests  are  shown  in  Fig.  27  and 
28.  These  plots  show  the  frequencies  of  the  two  pendulum  modes  and 
indicate  that  they  are  shorter  than  the  turning  period.  The  damping  of 
the  disturbances  is  so  low  that  one  might  assume  a slight  instability, 
if  eye  observation  had  not  indicated  that  the  disturbances  never  increased, 
no  matter  how  long  the  disturbances  were  observed. 

The  stability  being  so  small  one  would  expect  to  find  some  cases  of 
instability  if  the  weight  of  the  gondola,  the  weight  of  the  cable,  the 
drag  of  the  gondola,  and  the  turning  speed  were  varied.  For  all  these 
variations,  however,  no  change  of  stability  was  observed  within  the 
accuracy  of  the  measurements.  This  result  was  considered  a crucial 
and  sufficient  check  of  the  theory. 

The  rotary  motions  described  by  the  angles^q  and^_^  could  be 
observed  in  a small  number  of  cases,  (Fig.  26),  which  showed  that  these 
modes  have  a shorter  period  and  high  damping  which  causes  their  decay 
before  the  pendulum  modes  have  reached  their  first  maximum.  This  result 
justifies  the  separation  of  rotary  and  pendulum  modes  made  in  the 
theoretical  investigations. 

A detailed  evaluation  of  the  overhead  camera  pictures  was  considered 
not  justifiable.  Such  an  evaluation  would  have  required  a record  which 
includes  the  undisturbed  motion  and  a high  number  of  turns  after  the 
disturbance.  In  order  to  discern  between  the  different  turns  it  would 
have  been  necessary  to  install  a flicker  mechanism  in  the  circuit  of  the 
gondola  bulb  which  produces  different  signals  composed  of  dashes  and  dot3. 
This  installation  was  not  made  since  the  obtained  qualitative  checks  of 
the  theory  were  considered  sufficient. 
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Fig*  22  Record  Taken  by  the  Overhead  Camera 


Fig.  23  Record  Taken  by  the  Overhead  Camera  - 
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Fig.  2i|  Record  Taken  by  the  Overhead  Camera 


Fig*  25  Record  Taken  by  the  Floor  Camera 
WADC  TR  52-7S 
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Fig.  26  Record  Taken  by  the  Floor  Camera 
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Explanation  of  the  Observed  Instability 


Since  the  instability  observed  on  the  full  scale  equipment  could 
neither  be  reproduced  by  the  theory  for  a single  point  suspension  nor  by 
the  corresponding  model  tests,  the  reason  for  the  instability  must  be 
looked  for  in  the  difference  between  the  single  cable  and  the  double  cable 
suspension.  This  is  the  more  obvious  because  the  instability  at  50  mph 
did  not  develop  on  the  full  scale  equipment  any  more  after  the  double 
cable  suspension  had  been  replaced,  by  a single  cable  suspension.  From 
an  academic  standpoint,  an  exact  proof  of  this  statement  would  be  of 
interest. 

A numerical  investigation  of  this  instability,  however,  meets  the 
previously  mentioned  difficulties  caused  by  the  lack  of  the  necessary 
data  on  initial  twist  and  friction  between  the  sheathing  units.  Since 
the  Equipment  Laboratory  was  not  interested  in  such  an  investigation,  it 
was  not  undertaken. 


Conclusions  and  Recommendations 


Based  on  the  agreement  between  theory  and  tests,  the  following 
conclusions  are  drawn: 

1.  Regardless  what  the  weight  and  drag  of  the  gondola  and  of  the 

cable  might  be,  regardless  of  the  length  of  the  whirl  arm  and  of  the 
cable,  and  regardless  of  the  whirling  speed,  the  system  is  always  stable 
if  a one  point  suspension  is  used  and  if  Cp  - £ 0,  i.e.  if  the 

drag  is  higher  or  equal  to  the  slope  of  the  cross  force  coefficient 
plotted  versus  angle  of  yaw. 

2.  A violation  of  this  condition  mentioned  in  point  1 seems  to 
have  a minor  effect  on  the  stability. 

3.  The  observed  instability  was  very  likely  caused  by  irregularities 
in  the  mounting  of  cable  sheathing  and  by  mutual  interference  between 

the  aerodynamic  forces  on  the  two  cable  sheathings. 

4.  Because  of  the  observed  instability  and  because  of  the 
difficulties  to  be  expected  in  a theoretical  or  experimental  investi- 
gation of  the  double  point  suspension,  a single  point  suspension  is 
recommended. 

5.  From  a stability  standpoint  an  unsheathed  cable  of  circular 
cross-section  is  the  safest.  Since,  however,  the  high  drag  of  such  a 
cable  would  prevent  the  system  reaching  the  required  circumferential 
speed  of  500  mph  a sheathing  should  be  applied,  which  is  flutter  safe  and 
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has,  under  operational  conditions,  a torsional  frequency  high  enough 
not  to  couple  with  the  yawing  motion  of  the  gondola.  This  means  the 
mass  center  of  a unit  length  of  cable  and  sheathing  should  be  a few 
percent  of  the  sheathing  chord  in  front  of  the  center  of  pressure  and 
of  the  elastic  axis, 

6.  The  new  gondola  should  have  such  a form  that  its  slope  of  cross 
force  coefficient  versus  angle  of  yaw  is  as  small  as  possible.  This 
means  that  the  lateral  projections  of  all  surfaces,  which  produce  cross 
forces,  should  be  kept  to  a minimum.  If  there  would  be  no  drag  restrictions 
a sphere  would  be  the  ideal  form 
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TABLES 


Table  I 

Data  on  the  Full  Scale  Test  Equipment 
Used  in  Stability  Calculations 

Altitude  of  the  tower  120  ft 

Length  of  cables L = 116  ft 

Length  of  Protruding  Whirl  Am R * 56  ft 

Weight  of  Gondola  Including  Parachute  and  Dummy  . . . . mg  = 685  lbs 

Snpty  weight  of  Gondola  Including  Balancing  Lead  ....  405  lbs 

Load  of  Gondola 280  lbs 

Weight  per  ft  of  6 Stranded  Cables  Without  Sheathing.  . . = 5.17  lb/ft 

Weight  per  ft  of  6 Stranded  Cables  With  Sheathing.  . . . SPc.  = 6.27  lb/ft 

Length  of  Gondola 1=8. 5 ft 

Cross  Section  Area  of  Gondola . . ® = 3.14  ft^ 

Center  of  Gravity  of  Gondola  is  Located  at  41$  of  Length  of  Gondola 

Slope  of  Pitching  Moment  Coefficient  of  Gondola  ....  = 0.47 

Nondimensional  Pitching  Moment  Coefficient  *001075 

Drag  Coefficient  of  Gondola  Plus  Cable Cp  = 0.24 

Gondola  Mass  Factor  for  Loaded  Condition  ......  yU  = 335 

Mass  Density  of  Air  . . . 002378  lb  sec2 

ftA 

Maximum  circumferential  design  speed  . Vmax  = 500  mph 

Time  Factor » 2848.6 

vo 

Output  of  Driving  Motor  2300  HP 

Total  Gearing  Ratio  From  Motor  Shaft  to  Whirl  Am.  . . . 22:1 
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Table  II 


Numerical  Data  of  the  Model  Test  Equipment 


Length  Scale  Factor  . . 

Time  Scale  ....... 

Altitude  of  the  Tower  . 

Length  of  Cables  . . . . . 
Length  of  Protruding  Whirl  Arm 
Gondola  Weights  . ; . . . 


Weights  of  Cables 
Length  of  Gondola 


_/Y  - 1:20 

1:  -\J~2 O' 

6.0  ft 

L - 5.8  ft 

R « 2.8  ft 

Wr  - 17.9  gr 
W2  = 23.0  gr 
W3  - 38.8  gr 

WCl  - 17.1  gr  WC2-34.2gr 

5.6  inches 


Cross  Section  Area  of  Gondola 1.54  sq  in. 

Diameter  of  Drag  Disks D^«0 

D2  = 7/16" 
D3  - 7/8" 
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10  32800  1640  71.7  x 1C)6  415  x 1C*  22.7  12.8  10.0  8.5  7.5  .54  x 10"3  .06  x 10~2 

5 72400  3 510  468  x 106  1320  x 1CA  30.2  12.2  10.4  8.3  7.7  .13  x 10-8  .04  x 10“2 


